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ABSTRACT 

We present a mid-infrared investigation of the scaHng relations between supermas- 
sive black hole masses (Mbh) and the structural parameters of the host spheroids in 
local galaxies. This work is based on two-dimensional bulge-disk decompositions of 
Spitzer/IUAC 3.6 jiia images of 57 galaxies with Mbh estimates. We first verify the 
accuracy of our decomposition by examining the fundamental plane (FP) of spheroids 
at 3.6 /im. Our estimates of effective radii {Rq) and average surface brightnesses, com- 
bined with velocity dispersions from the literature, define a FP relation consistent 
with previous determinations but doubling the observed range in Rq. None of our 
galaxies is an outlier of the FP, demonstrating the accuracy of our bulge-disk decom- 
position which also allows us to independently identify pseudobulges in our sample. 
We calibrate M/L at 3.6 /im by using the tight M^yn-^bui relation (~ 0.1 dex intrinsic 
dispersion) and find that no color corrections are required to estimate the stellar mass. 
The 3.6 /im luminosity is thus the best tracer of stellar mass yet studied. We then ex- 
plore the connection between Mbh and bulge structural parameters (luminosity, mass, 
effective radius). We find tight correlations of Mbh with both 3.6 /im bulge luminosity 
and dynamical mass (MBn/Mdyn ~ 1/1000), with intrinsic dispersions of ~ 0.35 dex, 
similar to the Mbh— cr relation. Our results are consistent with previous determina- 
tions at shorter wavelengths. By using our calibrated M/L, we rescale MBH-^bui to 
obtain the Mbh-M^ relation, which can be used as the local reference for high-z stud- 
ies which probe the cosmic evolution of MBH-galaxy relations and where the stellar 
mass is inferred directly from luminosity measurements. The analysis of pseudobulges 
shows that 4 out of 9 lie on the scaling relations within the observed scatter, while 
those with small Mbh are significantly displaced. We explore the different origins for 
such behavior, while considering the possibility of nuclear morphological components 
not reproduced by our two-dimensional decomposition. 

Key words: black hole physics -galaxies: bulges - galaxies: fundamental parameters 
- galaxies: nuclei - galaxies: photometry - infrared: galaxies 



1 INTRODUCTION 

Supermassive black holes (BHs) are believed to dwell in 
almost all galaxy spheroids (hereafter bulges) as a result 
of past quasar activity (Soltan 1982, Marconi et al. 2004, 
Shankar 2009 and references therein) . A direct link between 
bulge formation and the growth of the central BHs has been 
inferred from tight relations between BH mass (Mbh) and 
the bulge structural parameters, such as velocity disper- 
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sion (cr, Ferrarese and Merrit 2000, Gebhardt et al. 2000), 
luminosity (Lbui, Kormendy & Richstone 1995, Marconi 
& Hunt 2003 - hereafter, MH03) and mass (Magorrian et 
al. 1998, MH03, Baring & Rix 2004, hereafter HR04). 
The underlying physics is encoded not only in the slopes of 
the correlations, but also in the intrinsic or cosmic scatter 
(rms), i.e. the scatter not accounted for by measurement 
uncertainties. Particular care has been taken with its mea- 
surement and, for instance, over the last 15 years the intrin- 
sic scatter of Mbh- Lhui decreased from an initial value of 
rms ^ 0.5 dex (in log Mbh, Kormendy & Richstone 1995) 
to - 0.3 - 0.4 dex (MH03, Giiltekin et al. 2009, Hu 2009 
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hereafter G09 and H09 respectively). Mbh-c" is considered 
one of the tightest relations {rms ^ 0.25 — 0.3 dex, Tremaine 
et al. 2002, T02), suggesting that bulge dynamics or mass 
rather than luminosity are driving the correlations. Indeed 
the MBH-Mdyn relation has a scatter comparable with Mbh- 
cr {r^ 0.3 dex, MH03, HR04, H09 and references therein). 
Nonetheless, slopes and rms of scaling relations are still a 
subject of debate; an accurate assessment of the actual rms 
is indeed crucial to draw conclusions on the physical origin 
of the correlations and turns out to be fundamental to both 
constrain BH-galaxy evolutionary models and compute the 
space density of BHs (e.g. G09, Hopkins et al. 2007). 

Since a successful theoretical model should be able 
to account for all empirical BH-bulge relations, discerning 
which are the truly fundamental ones, there has also been 
a growing interest in searching for all possible connections 
of the BH with other structural parameters. These include 
the effective radius Re (MH03) and its various combinations 
with a such as virial mass, gravitational binding energy 
of the host spheroid (e.g. MH03, Aller k Richstone 2007, 
H09, Feoli et al. 2010), thus leading to a possible fundamen- 
tal plane for BHs (Barway & Kembhavi 2007, Hopkins et 
al. 2007). 

To investigate the interplay between Mbh, bulge dy- 
namics and luminosity, it is crucial (z) to perform a 
proper photometric decomposition of the galaxy components 
(e.g. disk, bulge, bar etc.) and (ii) to obtain good estimators 
of stellar mass (M^), dynamical mass (Mdyn) and the galaxy 
mass-to- light ratio (M/L). 

An accurate bulge-disk decomposition can be performed 
by simultaneously fitting multiple bi-dimensional compo- 
nents to galaxy images. To compute the bulge mass, one can 
apply the virial theorem which turns out to be quite accurate 
(Cappellari et al. 2006) but requires a spectroscopical esti- 
mate of the velocity dispersion in addition to the bulge effec- 
tive radius Re- Alternatively, one can estimate by com- 
bining Lbui with a M/L ratio based on calibrated relations 
available in the literature (see, e.g., Bell & de Jong 2001, 
Beh et al. 2003, Cappellari et al. 2006). In this case, can 
be affected by dust extinction, depending on wavebands, 
and by the assumption of a constant M/L over the entire 
galaxy. In any case, it is believed that is more easily 
estimated than Mdyn especially at higher redshifts where it 
is the only host property not affected by strong biases (Mer- 
loni et al. 2010, Trakhtenbrot & Netzer 2010, Lamastra et 
al. 2010). 

Recently, attention has been drawn to the coexistence 
of pseudobulges and BHs (Graham 2008b, Greene, Ho & 
Barth 2008, H09, Nowak et al. 2010, Erwin 2010). Pseudob- 
ulges are bulges which are photometrically and morphologi- 
cally disk-like (possibly containing typical disk features such 
as bars, rings or ovals) and which present 'cool' kinematics, 
dominated by rotation (Kormendy & Kennicutt 2004). Ap- 
parently, pseudobulges follow their own relation with BHs, 
hosting less massive BHs than classical bulges (Greene, Ho, 
& Barth 2008, H09). Indeed, in these kind of structures, 
Mbh seems to better correlate with the small classical bulge 
component only (Nowak et al. 2010, Erwin 2010). However, 
these results are only based on a few pseudobulges 5) 
which have been analyzed as a separate class, and the poor 
statistics preclude a firm conclusion. 

In this paper we present a mid-infrared (MIR) view of 



the Mbh -bulge scaling relations based on a 2D photomet- 
ric decomposition of Spitzer/IRAC images at 3.6 //m. The 
superb performance of Spitzer/IRAC (Fazio et al. 2004) pro- 
vides images of unprecedented quality in the MIR for the 57 
galaxies analyzed here, and IRAC sensitivity permits the 
clear identification of morphological features. Therefore we 
identify classical bulges and pseudobulges with our analy- 
sis and compare their properties with Mbh for one of the 
largest samples ever considered. 3.6 fim observations are an 
ideal tracer of M^, and are less affected by extinction than 
shorter wavelengths. 
Our aims are to: 

- determine the bulge structural parameters with high accu- 
racy; 

- investigate the Mbh -bulge scaling relations and their in- 
trinsic scatter, taking into account the nature of the bulge 
(classical vs pseudobulges); 

- calibrate the M/L ratio thus supplying the reference Mbh- 
M^ in the local Universe for studies of the Mbh -galaxy scal- 
ing relations at higher z, based on M^ derived from luminos- 
ity measurements (by, e.g., Merloni et al. 2010, Trakhtenbrot 
& Netzer 2010). 

In Section 2 we describe sample selection, data reduc- 
tion and the grid-method adopted to decompose the IRAC 
3.6 fim images. We then present the fundamental plane (FP) 
for spheroids defined by the objects in our sample. Mbh- 
bulge scaling relations and the linear regression methods 
adopted to estimate slopes and intrinsic dispersion of the 
correlations are described in Section 3. Section 4 is dedi- 
cated to a discussion about (a) the utility of MnH-Mdyn as 
a benchmark to probe BH vs host galaxy co-evolution, and 
(b) pseudobulges location in the scaling relations. 

Throughout this paper we assume the standard cosmol- 
ogy with Ho = 70 km s"^ Mpc"\ Qm = 0.3, Qa = 0.7. 



2 DATA SELECTION AND ANALYSIS 

In this section we describe the criteria adopted to select our 
sources (Section 2.1), the techniques adopted to co-add the 
basic calibrated data (BCD) provided by version 18.7 of the 
Spitzer pipeline (Section 2.2), and the 2D fitting procedure 
followed to decompose IRAC 3.6 fim images (Section 2.3). 
Finally Section 2.4 contains the analysis of the FP at 3.6 
/im. 



2.1 Sample selection 

We consider the samples compiled by G09 and H09 for a to- 
tal of 64 galaxies. These authors have selected only galaxies 
with reliable Mbh estimates (i.e. by stellar or gas dynamics 
and masers). This is one of the largest samples used so far 
to investigate BH-galaxy scaling relations which also allows 
a consistent comparison of our results with previous works. 

G09 and Batcheldor (2010) have recently shown that 
discarding sources whose sphere of influence is not spatially 
resolved can bias or even simulate the observed linear re- 
lations. Therefore, we do not apply any cut based on that 
criterion. To increase the number of objects with low Mbh, 
we add to the sample NGC3368 and NGC3489 whose Mbh 
is estimated in Nowak et al. (2010). These objects have also 
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OBSERVATIONS 



Source Name 


PID 


Date 


Frame Time 


Frames 


Position 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Circinus 


40936 


2007-09-09 


12 


1 


4 


IC1459 


20371 


2005-11-28 


2 


1 


22 


IC2560 


40936 


2007-07-05 


12 


1 


4 


IC4296 


20371 


2006-02-13 


2 


1 


22 


NGC221 


00069 


2004-07-19 


12 


1 


5 


NGC524 


50630 


2008-09-18 


30 


1 


5 


NGC821 


20371 


2005-08-21 


2 


1 


22 


NGC1023 


00069 


2004-02-11 


30 


1 


5 


NGC1068 


00032 


2004-01-16 


12 


2 


5 


NGC1300 


61065 


2009-09-06 


30 


1 


11 


NGC1316 


00159 


2004-02-23 


30 


1 


9 


NGC2549 


50630 


2008-12-23 


30 


1 


5 


NGC2748 


61063 


2009-12-02 


30 


1 


4 


NGC2778 


30318 


2006-11-24 


100 


1 


5 


NGC2787 


03674 


2004-10-30 


30 


1 


3 


NGG2974 


30318 


2006-12-28 


100 


1 


5 


NGC3031 


00159 


2004-05-01 


30 


1 


12 


NGC3079 


00059 


2004-04-05 


12 


1 


4 


NGC3115 


00069 


2004-04-29 


12 


1 


4 


NGC3227 


03269 


2004-12-21 


12 


1 


1 


NGC3245 


03674 


2004-11-28 


30 


1 


3 


NGC3368 


00069 


2004-05-19 


30 


1 


5 


NGC3377 


00069 


2004-05-27 


30 


1 


5 


NGC3379 


00069 


2004-12-15 


12 


1 


5 


NGC3384 


30318 


2006-12-27 


100 


1 


5 


NGG3414 


50630 


2009-01-29 


30 


1 


5 


NGC3489 


00069 


2004-05-19 


12 


1 


5 


NGC3585 


30318 


2006-07-09 


100 


1 


5 


NGC3607 


00069 


2004-05-19 


12 


1 


5 


NGC3608 


30218 


2006-12-27 


100 


1 


5 


NGC3998 


00069 


2004-04-21 


12 


1 


5 


NGC4026 


30318 


2006-12-26 


100 


1 


5 


NGC4151 


03269 


2004-12-17 


12 


1 


5 


NGC4258 


20801 


2005-12-25 


30 


2 


7 


NGC4261 


00069 


2004-05-27 


12 


1 


5 


NGC4374 


00069 


2004-05-27 


12 


1 


5 


NGC4459 


03649 


2005-01-22 


12 


1 


5 


NGC4473 


03649 


2005-01-22 


12 


1 


5 


M87" 


03228 


2005-06-11 


30 


1 


5 




03228 


2005-06-11 


30 


1 


5 


NGC4486A 


03228 


2005-06-11 


30 


1 


5 


NGC4552 


00159 


2004-05-27 


30 


1 


7 


NGC4564 


20371 


2006-02-09 


2 


1 


22 


NGC4594 


00159 


2004-06-10 


30 


1 


6 


NGG4596 


03674 


2005-06-10 


30 


1 


3 


NGC4621 


03649 


2005-06-10 


12 


1 


5 


NGC4649 


00069 


2004-06-10 


12 


1 


5 


NGC4697 


03403 


2005-06-17 


30 


1 


5 


NGC5077 


00069 


2005-05-11 


12 


1 


5 


CenA 


00101 


2004-02-11 


12 


1 


5 


NGC5576 


03403 


2005-07-15 


30 


1 


5 


NGC5813 


00069 


2004-02-17 


12 


1 


5 


NGC5845 


20371 


2005-08-23 


2 


1 


22 


NGC5846 


00069 


2004-03-09 


12 


1 


5 


NGG6251 


02418 


2004-12-16 


30 


2 


6 


NGC7052 


30877 


2006-11-26 


30 


1 


5 


NGC7457 


30318 


2006-07-12 


100 


1 


5 


NGC7582 


03269 


2004-11-27 


12 


1 


1 



Table 1. Observations Log. Columns: (1) Source name. " mean of 2 observations (see Section 2.2 for details). (2) Proposal identification 
number (PID). (3) Observing date. (4) Exposure time for each frame. (5) Numbers of frames. (6) Number of positions to realize the 
source maps. 
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a pseudobulge and their inclusion allows us to better investi- 
gate at 3.6 /im the discrepancies between pseudo and classi- 
cal bulge properties found by H09 and Nowak et al. (2010). 

We thus obtain a sample of 66 galaxies, 7 of which 
(CygnusA, NGC1399, NGC3393, NGC4291, NGC4342, 
NGC4742, NGC5252) have no 3.6 /im public Spitzer/IRAC 
observations. In addition, we exclude M31 because the avail- 
able maps do not cover the entire galaxy and we exclude the 
Milky Way as there are no secure measurements of its 3.6 
/im bulge properties. The final sample thus consists of 57 
galaxies. Tables [l] and |2] provide the basic parameters of 
the observations together with distances, Mbh and a mea- 
surement Q To identify pseudobulges in 3.6 fim images we 
adopt a criterion based on the Sersic index and examine a 
posteriori the sources location in the 3.6 /im FP of elliptical 
galaxies. 

2.2 IRAC images 

Given the wavelength dependence of the emission from stel- 
lar photospheres, the 3.6 /xm IRAC band is the best suited 
for our study. In fact, the STARBURST99 models (Leitherer 
et al. 1999) show that only ^ 10% and even less of the 
emission in the 4.5 fim band is stellar, and that at longer 
wavelengths the fraction is even lower (Helou et al. 2004). 
Hence we used only channel 1 (centered at 3.55 /xm) basic 
calibrated data (BCD). 

Individual BCD frames were processed with version 18.7 
of the SSC pipeline. The image mosaicking and source ex- 
traction package (MOPEX, Makovoz & Marleau 2005) was 
used to co-add BCD frames for each source. Bad pixels, 
i.e. pixels masked in the Data Collection Event (DCE) sta- 
tus files and in the static masks (pmasks), were ignored. 
Additional inconsistent pixels were removed by means of the 
MOPEX outlier rejection algorithms. We relied on the dual 
outlier technique, together with the multi-frame algorithm. 
The frames were corrected for geometrical distortion and 
projected onto a North-East coordinate system with pixel 
sizes of 1.20 arcsec, equivalent to the original pixels. Mosaics 
for data images were realized with standard linear interpola- 
tion. The same was done for the uncertainty images, i.e. the 
maps for the standard deviations of the data frames. The 
signal-to-noise ratio in our post-pipeline MOPEX mosaics 
are comparable to or better than those in the SSC products. 
M87 has two maps with different spatial extension, thus we 
average them with a weighted mean where the weights are 
given by uncertainty images. 

2.3 Image analysis 

We use the public program GALFIT V.3 (Peng, et al. 2002, 
2007) to perform a two-dimensional decomposition of 3.6 
/xm images in bulge/disk components. Besides the standard 
inputs (e.g. data, point spread functior[^(PSF) images etc.), 

^ We assume that the effective velocity dispersion ae is consis- 
tent, within the errors, with the central velocity dispersion ac and 
we refer to it as a. This choice is supported by G09, whose Fig. 2 
shows that no bias is implied by this choice. 

^ We use the instrumental PSF image provided by the SSC, after 
checking that its full width at half maximum (FWHM) is consis- 
tent with what observed for the foreground stars in our images. 



GALFIT requires a standard-deviation image, used to give 
relative weights to the pixels during the fit, and a bad pixel 
mask. We therefore use the uncertainty data obtained from 
MOPEX as sigma images and construct a bad pixel frame 
masking out foreground stars, background galaxies and pos- 
sible irregularly shaped regions such as dust lanes across the 
galaxy. 

We choose the number and kind of model components 
on the basis of the Hubble morphological types (listed in 
Table [T]) and after a visual inspection of the images. Thus, 
for elliptical galaxies we use a pure Sersic profile, while for 
lenticular (SO) and spiral galaxies we add an exponential 
disk. In the case of barred galaxies (SB), we consider an 
additional Gaussian component, equivalent to a Sersic pro- 
file with index n — 0.5. In some cases, even if the source 
is classified as a barred galaxy, the bar cannot be identified 
in the MIR (e.g. NGC2778). In these cases we do not add 
any Gaussian component to the model. In the case of active 
galaxies, we also include a nuclear point source to account 
for the emission of the active galactic nucleus (AGN). 

The contribution of thermal dust to the observed emis- 
sion at 3.6 /im can be powered by an active galactic nucleus 
(AGN) and/or star formation. The former can be accounted 
for by including a nuclear PSF in the model, while the lat- 
ter can be masked out. Indeed intense star formation occurs 
in well- localized knots of emission, and does not affect the 
entire galaxy disk. Overall we estimate that hot dust can 
contribute no more than a few percent to the bulge (or disk) 
total luminosity and, given its localized nature, should not 
significantly affect the structural parameters derived from 
our global two-dimensional fits. 

We fix the background in the fits, estimating it as the 
mean surface brightness (with the relative standard devia- 
tion) over an annular region surrounding the galaxy between 
2-3 times the optical radius. Foreground sources such as stars 
or galaxies are not considered in the background calculation 
by means of a 2.5 sigma rejection criterion. 

Some aspects of the photometric decomposition are 
potentially problematic. Due to the functional form of 
the Sersic profile, the effective radius Re is coupled to 
the Sersic index n. To control this degeneracy we fol- 
low Hunt, Pierini & Giovanardi (2004, hereafter HPG04) 
and fit every galaxy with Sersic indexes n fixed to values 
n = 1, 1.5, 2, 2.5, 3, 3.5, 4, 5, 6, 7, which reasonably span all 
galaxy morphologies. We then assume 0.5 as the error on 
the Sersic index retrieved from the best fit (see below for 
the criteria defining the best model) . 

In addition, background measurements can be affected 
by galaxy extended emission; indeed, in several cases the 
faint tail of source emission covers most of the field-of-view 
(NGC3115, NGC4621, NGC4374, NCG5846, NGC7052 are 
some extreme examples) . Thus for each galaxy, we fix three 
different background values: s^s — Gs and s + cts, where s and 
Gs are the sky flux and its standard deviation estimated as 
described above. This allows to account for possible large 
variations of the background. 

Once the number of components is defined, we run 30 
different models; for each of the 10 fixed n values we use the 
3 different background estimates. The free parameters are 
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SAMPLE PROPERTIES. 



Source Name 


Type 


Distani 


(1) 


(2) 


(3) 


Circinus 


Sb 


4.0 


IC1459 


E4 


30.9 


IC2560 


SBb 


40.7 


IC4296 


E 


50.8 


NGC221 


E2 


0.86 


NGC524 


SO 


33.6 


NGC821 


E4 


25.5 


NGC1023 


SBO 


12.1 


NGC1068 


Sb 


15.4 


NGC1300 


SB(rs)bc 


20.1 


NGC1316 


SBO 


19 


NGC2549 


SO 


12.3 


NGC2748 


Sc 


24.9 


NGC2778 


E2 


24.2 


NGC2787 


SBO 


7.9 


NGC2974 


E4 


21.5 


NGC3031 


Sb 


4.1 


NGC3079 


SBcd 


15.9 


NGC3115 


SO 


10.2 


NGC3227 


SBa 


17.0 


NGC3245 


SO 


22.1 


NGC3368 


SABa 


10.4 


NGC3377 


E6 


11.7 


NGC3379 


EO 


11.7 


NGC3384 


SBO 


11.7 


NGC3414 


SO 


25.2 


NGC3489 


SBO 


12.1 


NGC3585 


SO 


21.2 


NGC3607 


El 


19.9 


NGC3608 


El 


23 


NGC3998 


SO 


14.9 


NGC4026 


SO 


15.6 


NGC4151 


Sa 


20.0 


NGC4258 


SABbc 


7.2 


NGC4261 


E2 


33.4 


NGC4374 


El 


17 


NGC4459 


E2 


17 


NGC4473 


E5 


17 


M87 


El 


17 


NGC4486A 


E2 


17 


NGC4552 


E 


15.3 


NGC4564 


SO 


17.0 


NGC4594 


Sa 


10.3 


NGC4596 


SBO 


18.0 


NGC4621 


E5 


18.3 


NGC4649 


E2 


16.5 


NGC4697 


E6 


12.4 


NGC5077 


E3 


44.9 


CenA 


SO/E 


4.4 


NGC5576 


E3 


27.1 


NGC5813 


El 


32.2 


NGC5845 


E3 


28.7 


NGC5846 


EO 


24.9 


NGC6251 


El 


106 


NGC7052 


E3 


70.9 


NGC7457 


SO 


14.0 


NGC7582 


SBab 


22.0 



(4) 



(5) 



logLy 

(6) 



logLx 
(7) 



7+0.004 
-0.003 
8+11 
-12.0 
044+0-°^^ 



0.017^ 
28^ 



13.5 



-0.022 
+ 2 



031+° °°^ 
^+0.60 
-0.37 



8.32^ 



0.42 



+ 0.28 
-0.08 



-0.003 
+0.69 
-0.35 

^••^^-0.27 
^+0.10 
-0.12 
7+0.38 
-0.38 

^ "-0.10 

^^^-0.05 
1 7+0.3 
-0.3 

o-8ol°"ii 

0.025lg-- 

R + 5-4 
-2.9 
r,+0.37 
-0.08 
^ + 0.5 

0.075+r.^'. 



0.71^ 



0.14^ 
0.47^ 



9.6^ 
0.20^ 
2.2^ 



1.1 
1.2 
0.18 
2.51 



0.015 
+ 1.1 
-0.1 



9 + 0. J 

-0.58 
e + 0.01 
-0.03 
-, +0.30 
-0.31 
06+°-°°^2 
-0.0012 



.+1.5 



1.2^ 



3.4^^ o 
9+0^^4° 
-0.41 

2 1 + 1-1 
-0.7 
2 4+2.1 

1+0.7 
-0.4 
^ + 0.07 
-0.07 
0.378+0-00; 



2-i: 

0.65^ 



5.5 
15 
0.74 



0.001 
+ 1.1 
1.2 
+ 11 



i+0.14 
-0.14 

1 q+0.5 

l-'^-0.94 

36 + 1° 
'^'^-lO 
.0+0.05 
^•l'^-0.04 
5 + °-^ 

0.691-- 



0.84' 



4.0 
+0.36 
0.25 

^•^-0.4 
21+! 

2.0^ 

8.0^ 

3.0^ 

lA 



n+0.2 
-0.2 

n+5.0 

-3.3 

n+0-4 

-0.2 
« + 0.3 
-0.4 
7 1+0.9 
-0.6 
Q + 0.5 
-1.7 

6.01-0 

041+°-°12 
-0.017 

p; + 0.26 

-0.19 



2.9^ 



0.55^ 



158±18 
340±17 

144±7 
226±10 

75 ±3 
235 + 12 
209 + 10 
205 + 10 

151 + 7 
218+10 
226+11 

145+7 

115+5 

175+8 

189+9 
227+11 

143 + 7 

146 + 7 
230+11 

133+6 
205+10 
98.5+5 

145+7 
206+10 

143+7 
205+10 

91 +5 
213 + 10 
229 + 11 

182+9 
305+15 

180+9 

156+8 
115+10 
315+15 
296 + 14 

167+8 

190 + 9 
375 + 18 

111 + 5 
252+12 

162+8 
240+12 

136+6 
225+11 
385+19 

177+8 
222 + 11 

150 + 7 

183+9 
230+11 
234+11 
237+12 
290+14 
266+13 

67 +3 
156 + 19 



O+0.06 
°-0.06 



R+0.02 
°-0.02 



+ 0.05 
0.06 
+ 0.10 
0.13 
-1^ 4. + 0.07 
10-^-0.09 



10.42 
10.17 



4.+0.05 
-0.06 



4+0.02 
-0.02 



9 97+0.04 

-0.04 

10 37+0-01 

9 90+°-°^ 
^•^^-0.10 



^ + 0.07 
^-0.09 



+ 0.04 
r+0.04 
^-0.04 



R+0.07 
°-0.09 



11.02"^ 



-0.02 
10 91 + 0-02 



11.1 

9.41 



-0.02 
+ 0.02 
-0.02 
+ 0.02 
0.02 



9.77^ 



7+0.11 
-0.15 
10 9+0-06 
1^-^-0.06 



10.99 

10.44 

10.74 

10. 

10.43 



+ 0.02 
0.02 
+ 0.04 
0.05 
+ 0.05 
0.06 
+ 0.03 
-0.03 
+ 0.05 
0.06 



.+0.05 
-0.06 



9 42+0-04 
-0.05 



10.03 
11.59 
10.46 
12.36 

8.83 
11.46 

10.8 
10.54 
10.82 

11.2 
10.18 



9.86 
10.95 
10.44 
10.26 
10.46 
9.92 
10.69 

10.22 
10.96 
10.44 
10.71 

11.18 
11.09 
10.81 
10.63 
10.57 
10.27 
9.93 
11.42 

10.58 
10.92 
11.46 
10.19 
11.04 
10.35 

10.44 
10.77 
11.46 
10.52 
11.35 
10.49 

10.5 
11.06 
10.62 
11.33 
11.81 
11.39 

9.69 



Table 2. Columns: (1) galaxy name. (2) Hubble type taken from the Hyper-Leda catalogue. (3) Distance in Mpc. (4) Mbh in lO^M©, 
see G09, H09 and Nowak et al. (2010) for references on BH mass measurements. (5) Velocity dispersion in km/s. (6)-(7) logarithmic V- 
and K-band luminosities in solar units. The la error of Lx is 10% (H09). All values in columns (4)-(7) are from G09 and H09 if not 
included in V-band sample, and are scaled from the original publications to our preferred distances (column 3) assuming the standard 
cosmology. 
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DECOMPOSITION PARAMETERS 



Source Name 


77, 


'"'OuLge 


Re 


(b/a)h 


Co 










FWHM 


(h / n?\u„^ 


777, A AT 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Circinus 


2 


5.48± 0.04 


0.21±0.01 


0.68 


0.2 


4.4 


1.58 


0.39 


- 


- 


- 


- 


IC1459 


6 


6.15± 0.3 


9.15±3.71 


0.76 


-0.14 


- 


- 


- 


- 


- 


- 


12 


IC2560 


2 


9.1 ± 0.7 


5.43±2.84 


0.37 


0.62 


9.0 


6.18 


0.57 


- 


- 


- 


10.8 


IC4296 


4 


7.09± 0.5 


8.28±5.91 


0.92 


0.66 


- 


- 


- 


- 


- 


- 


11.6 


NGC221 


4 


5.2 ± 0.63 


0.12±0.07 


0.7 


0.08 


6.5 


0.09 


0.89 


- 


- 


- 


- 


NGC524 


3 


7.01± 0.19 


4.37±2.54 


0.95 


-0.07 


8.7 


4.65 


- 


- 


- 


- 


- 


NGC821 


7 


6.9 ± 0.26 


7.86±4.15 


0.62 


-0.19 


- 


- 


- 


- 


- 


- 


- 


NGC1023 


3 


6.7 di 0.18 


1.41±0.80 


0.61 


-0.55 


6.67 


48.51 


0.3 


9.01 


1.13 


0.55 


- 


NGC1068 


1 


6.15ib 0.05 


0.73ib0.02 


0.73 


0.19 


6.1 


1.76 


0.81 


- 


- 


- 


7.37 


NGC1300 


3 


8.40± 0.33 


8.32 ±1.30 


0.29 


0.0 


7.56 


7.21 


0.74 


- 


- 


- 


- 


NGC1316 


5 


4.7 ± 0.8 


8.57±0.55 


0.71 


-0.15 


6.91 


11.52 


0.63 


- 


- 


- 


11.13 


NGC2549 


7 


8.32± 0.13 


0.69 ±0.15 


0.5 


-0.42 


9 


1.81 


0.24 


- 


- 


- 


- 


NGC2748 


4 


9.64± 0.09 


1.87 ±0.54 


0.95 


-0.31 


8.72 


1.9 


0.24 


- 


- 


- 


- 


NGC2778 


2.5 


10.44±0.15 


0.29±0.05 


0.73 


0.09 


9.74 


1.22 


0.79 


- 


- 


- 


- 


NGC2787 


3 


7.57± 0.57 


0.60±0.19 


0.62 


0.0 


8.07 


1.12 


0.52 


11.5 


0.17 


0.11 


12.3 


NGC2974 


3 


7.28± 0.14 


2.83±0.61 


0.67 


-0.06 


- 


- 


- 


- 


- 


- 


11.8 


NGC3031 


3 


4 ± 0.31 


2.53±0.87 


0.62 


-0.03 


4.51 


3.35 


0.49 


- 


- 


- 


10.2 


NGC3079 


2 


6.87ib 0.08 


5.71±0.37 


0.17 


0.33 


9.48 


5.31 


0.86 


10 


0.27 


0.13 


9.64 


NGC3115 


3 


5.99± 0.1 


1.35±0.15 


0.33 


-0.62 


6.69 


3.4 


0.41 


- 




- 


- 


NGC3227 


4 


7.27± 0.5 


6.83±2.27 


0.49 


0.18 


9.44 


2.36 


0.44 


- 




- 


- 


NGC3245 


2.5 


8.73± 0.2 


0.49±0.04 


0.58 


-0.13 


8.14 


2.3 


0.52 


- 




- 


- 


NGC3368 


1 


6.57± 0.1 


2.3 ±0.12 


0.63 


-0.13 


6.83 


6.17 


0.58 


7.25 


0.7 


0.55 


- 


NGC3377 


6 


6.93± 0.18 


3.13±0.88 


0.55 


0.03 


- 


- 


- 


- 


- 


- 


- 


NGC3379 


5 


5.85± 0.16 


2.61±0.64 


0.88 


0.03 


- 


- 


- 


- 


- 


- 


- 


NGC3384 


2.5 


8.11± 0.07 


0.25±0.01 


0.76 


-0.04 


7.06 


2.5 


0.49 


8.73 


1.13 


0.79 


- 


NGG3414 


5 


7.87± 0.5 


2.67 ±1.29 


0.75 


-0.44 


9.28 


3.04 


0.87 


12.4 


0.29 


0.36 


- 


NGC3489 


1.5 


8.3 lb 0.45 


0.27±0.21 


0.82 


0.15 


7.7 


1.11 


0.54 


- 


- 


- 


- 


NGC3585 


2.5 


7.21± 0.63 


1.59±1.97 


0.55 


-0.28 


7.08 


5.25 


0.61 


- 


- 


- 


- 


NGC3607 


5 


6.44± 0.25 


4.3 ±2.02 


0.83 


0.0 


- 


- 


- 


- 


- 


- 


12.9 


NGC3608 


6 


7.19± 0.23 


6.29±2.51 


0.79 


0.12 


- 


- 


- 


- 


- 


- 


- 


NGC3998 


1.5 


8.14± 0.05 


0.34±0.36 


0.82 


0.07 


7.75 


1.51 


0.79 


- 


- 


- 


10.3 


NGC4026 


3.5 


8.08± 0.18 


0.86±0.30 


0.60 


0.0 


8.36 


2.34 


0.16 


- 


- 


- 


- 


NGC4151 


3.5 


8.13± 0.41 


0.52±0.42 


0.87 


0.07 


8.17 


2.99 


0.61 


- 


- 


- 


14.6 


NGC4258 


2 


5.27di 0.08 


3.9 ±0.42 


0.46 


0.18 


8.73 


6.84 


0.66 


- 


- 


- 


9.6 


NGC4261 


4 


7.23d= 0.15 


3.66±2.80 


0.80 


0.0 


8.81 


6.65 


0.99 


- 


- 


- 


13.09 


NGC4374 


7 


5.06ib 0.04 


8.73±1.00 


0.90 


0.0 


- 


- 


- 


- 


- 


- 


11.7 


NGC4459 


2.5 


7.7 ± 0.1 


0.85±0.10 


0.83 


0.02 


7.7 


2.66 


0.76 


- 


- 


- 


- 


NGC4473 


7 


6.59± 0.05 


4.06±0.87 


0.56 


-0.07 


- 


- 


- 


- 


- 


- 


- 


M87 


4 


5.10± 0.05 


8.20±0.60 


0.82 


0.0 


- 


- 


- 


- 


- 


- 


12.26 


NGC4486A 


2.5 


9.04± 0.08 


0.59±0.07 


0.8 


0.03 


- 


- 


- 


- 


- 


- 


- 


NGC4552 


4 


6.65± 0.08 


1.8 ±0.81 


0.92 


-0.05 


7.78 


7.17 


0.8 


- 


- 


- 


11.5 


NGC4564 


7 


7.76± 0.13 


2.06±0.03 


0.62 


0.33 


9.46 


1.39 


0.2 


- 


- 


- 


- 


NGC4594 


1.5 


5.51dz 0.14 


3.3 ±0.35 


0.26 


-0.93 


5.79 


2.74 


0.61 


6.32 


7.79 


0.71 


9.64 


NGC4596 


3 


7.6 ih 0.05 


2.44±0.12 


0.53 


-0.23 


7.66 


4.91 


0.83 


11.76 


0.05 


0.99 


- 


NGG4621 


5 


6.22± 0.1 


5.46±1.09 


0.63 


-0.19 


- 


- 


- 


- 


- 


- 


- 


NGC4649 


3 


5.55di 0.10 


3.77±0.34 


0.82 


0.0 


8.3 


3.3 


0.93 










NGC4697 


5 


5.69± 0.09 


6.04±1.39 


0.65 


-0.1 
















NGC5077 


6 


7.61± 0.26 


6.35±2.34 


0.71 


0.11 














15.53 


CenA 


3.5 


3.49± 0.27 


2.21±1.2 


0.83 


0.0 


4.16 


3.92 


5.20 


0.50 






10.2 


NGC5576 


7 


7.26± 0.15 


4.51±1.29 


0.71 


0.03 
















NGC5813 


6 


6.52± 0.18 


17.4±3.80 


0.78 


0.1 
















NGC5845 


3 


8.96± 0.02 


0.51±0.02 


0.67 


-0.06 
















NGC5846 


3 


6.70± 0.06 


4.40±0.90 


1 


0.0 
















NGC6251 


7 


8.2 ± 0.19 


21.8±7.71 


0.84 


-0.09 














12.8 


NGG7052 


5 


7.96di 0.16 


13.5±3.81 


0.49 


0.46 
















NGC7457 


7 


9.54± 0.90 


0.94±0.81 


0.52 


0.19 


8.37 


1.61 


0.56 










NGG7582 


4 


7.33± 0.44 


9.22±4.39 


0.5 


0.23 


8.14 


3.18 


0.22 








8.62 



Table 3. Columns: (1) galaxy name. (2)-(6) bulge structural parameters: Sersic index, 3.6 /im magnitude, effective radius in kpc, axis 
ratio, diskyness-boxiness (Co < for disky ellipsoids, Co > for the boxy ones). (7)-(9) disk properties: 3.6 /im magnitude, scale length 
in kpc, axis ratio. (10-12) bar properties: 3.6 /xm magnitude, FWHM in kpc, axis ratio. (13) AGN 3.6 fim magnitude. 
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the bulge and disk brightnes^ their scale lengths, elliptici- 
ties, position angles and the bulge diskyness/boxyness. Ad- 
ditional free parameters can be the bar brightness, FWHM, 
ellipticity and position angle and/or the AGN point-source 
3.6 /im magnitude. As starting guesses we use the axis ra- 
tio, the major axis position angle, provided by the Hyper- 
Ledsj^ database. The initial disk and bulge scale lengths {Rq 
and Re respectively), are proportional to the optical radius 
{Ropt listed in HyperLeda): Rq ^ 0.25i?opt; the effective ra- 
dius follows the bulge shape according to the Sersic index 
Re ~ f^Ropt for n < 4 and Re ^ OARopt for n ^ 4 (Giova- 
nardi & Hunt 1988, Moriondo, Giovanardi & Hunt 1998). 

As shown by HPG04, the bulge component influences 
the value only in the inner region of the galaxy; indeed, 
the reduced depends on n only at small distances (radii) 
from the center, while the large area of the disk strongly di- 
lutes the influence of the bulge on the global • Appendix A 
(available online) shows examples of the trend reduced 
with radii for our galaxies outlining that fits with acceptable 
X^ values can differ significantly in the inner region where 
the bulge is important. Therefore, following HPG04, we se- 
lected the fits with the best n and sky values considering 
the lowest reduced within half optical radius Ropt,o.5- 
The fits are then considered reliable on the basis of circum- 
nuclear flux conservation; the ratio between the inner counts 
in the residual and source images must lower than 0.05. This 
means that less than 5% of the central flux can be lost or 
introduced by the model. As described in the following sec- 
tion, the recovery of a tight FP at 3.6 /im shows the accuracy 
of our analysis in estimating bulge structural parameters. 

The absolute errors on the free parameters are derived 
as the variations between the best fit values and the ones ob- 
tained from the model having the closest s, s -\- as oi s — as 
value and the lowest possible reduced within Ropt, 0.5 ac- 
cording to the criterion described above. With this choice we 
overestimate the statistical fit errors, but carefully constrain 
the effects due to an uncertain estimate of the background. 
The fit parameters and their uncertainties are listed in Ta- 
ble [3] See also Appendix A (available online) for further 
details on the error estimates. 

2.4 The fundamental plane at 3.6 fim and the 
identification of pseudobulges 

Early type galaxies and the bulges of disk galaxies follow 
the well known relations connecting their effective radius 
Re, mean surface brightness < le > within Re, and veloc- 
ity dispersion a (e.g. Djorgovsky & Davis 1987, Prugniel & 
Simien 1996, Jorgensen, Franx & Kjaergaard 1996, Burstein 
et al. 1997, Cappellari et al.2006). In particular, Jun & Im 
(2008, hereafter JI08) studied the FP at 3.6 /im considering 
a sample of elliptical galaxies spanning an order of magni- 
tude in effective radii. Here, to analyze the bulge properties 
and verify the reliability of our analysis, we adopt the rela- 
tion found by JI08: 

log Re = 1.55 X logcr - 0.89 x log < h > -9.89, (1) 



where a in km/s. Re is in kpc and < /e > in mag arcsec"^. 

Figure [1] shows the 3.6 /xm FP for all our decomposed 
bulges obtained by using equation ([T]) (dashed black line). 
From Fig. [l] we can draw two main conclusions: (i) our 
bulges follow the 3.6 /xm FP of elliptical galaxies with low 
dispersion, and (ii) pseudobulges (open red squares) lie on 
the FP within the observed dispersion. Indeed, (i) if we 
follow JI08 and apply a maximum likelihood method with 
a bootstrap resampling to account for errors, the intrinsic 
dispersion is: 0.20 ± 0.04 dex for the entire sample, and 
0.11 di 0.04 dex for no n- active galaxies. The last result is 
consistent within the errors with JI08 (rms = 0.09 =b 0.04, 
private communication). Thus the presence of an AGN can 
slightly worsen the 2D decomposition, introducing some de- 
generacy with the bulge component. Nevertheless, the low 
dispersion of the 3.6 /xm FP for our bulges and the absence 
of strong outliers suggests that our 2D decomposition is reli- 
able. A tight FP and its extension to reliably compact bulges 
then permit us to calibrate M/L (Fig.[4]and see Section 3.2). 

A detailed analysis of the FP is beyond the scope of this 
paper. However, it should be noted that our sample spans an 
order of magnitude larger range in galaxy sizes, extending 
the FP to smaller galaxy radii than in JI08. This could be 
the origin of the possible small disagreement between the 
FP defined by our sources and the determination by JI08 
(see Fig. [1]). 

Classical bulges and pseudobulges differ in morphol- 
ogy and dynamics: the first behave like early- type galax- 
ies, scaled to small sizes, surrounded by a disk component. 
Classical bulges thus are either prolate or oblate spheroids 
and are kinematicall}j^ hot, (tending to be pressure domi- 
nated with a stellar V/a < 1). On the contrary pseudob- 
ulges are disk-like components with V/a > 1 cold kinemat- 
ics dominated by rotation and possibly containing disk fea- 
tures such as nuclear bars and rings as well as young stellar 
populations (Andredakis, Peletier & Balcells 1995, Fisher & 
Drory 2010, hereafter FDIO, and references therein). A cor- 
relation between the bulge morphology and the Sersic index 
has been proposed by considering the morphological prop- 
erties of pseudobulges as shown in FDIO. We thus follow the 
FDIO Sersic index criterion to classify those bulges in spi- 
rals and lenticulars with n < 2 as pseudobulges. By means 
of the Sersic index, we identify nine disk galaxies as host- 
ing pseudobulges: Circinus, IC2560, NGC1068, NGC3079, 
NGC3368, NGC3489, NGC3998, NGC4258, NGC4594. The 
pseudobulges (ii) which we have identified are marked with 
open red squares in Fig. [l] and are not outliers of the 3.6 
/im FP. This is not surprising since the physical meaning 
of the fundamental plane is that galaxies are virialized sys- 
tems (Burstein et al. 1997) and this applies also to rotation- 
ally supported pseudobulges. Our sample partly overlaps the 
FDIO one, and for the 8 sources in common the classical ver- 
sus pseudo classification of bulges is in perfect agreement. 



^ To compute the magnitudes at 3.6 fim we adopt a zero point 
of 17.25 in Vega magnitudes, according to the IRAC photometric 
system (Reach et al. 2005). 
^ http://leda.univ-lyonl.fr/ 



^ The bulge kinematics are defined according to Erwin (2010) 
where V is the stellar rotational velocity deprojected to its in- 
plane value and a is the stellar velocity dispersion. 
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Figure 1. MIR fundamental plane for our decomposed bulges. The 3.6 /im FP coordinate axes are given by equation ([T| (JI08), 
represented by the black dashed line. The bulge parameters are the effective radius Re in kpc (Tab.jsJ, the velocity dispersion a in km/s 
(Tab.jSJ and the mean surface brightness within Re, < lie >= ''^bui^e + 2-5 log[27r(6/a)5i?g], where {b/a)t is the bulge axis ratio. < /ie > is 
expressed in mag/arcsec^ (for a consistent comparison with JI08, Vega magnitudes in Tabjsjare converted to the AB photometric system 
as m(AB)=m(Vega)+2.79). The color codes reflect the adopted fltting model, i.e. black points correspond to pure Sersic proflles, blue 
symbols refer to Sersic and an exponential disk (open circles), possibly with an additional Gaussian component (asterisks). The empty 
green triangles correspond to active galaxies, whose AGN emission is reproduced adding a point-like source to the models listed above. 
Open red squares highlight bulges with n ^ 2 which we thus classify as pseudobulges according to the FDIO criterion; no additional 
components are added to model these galaxies. The gray-shaded region highlights the range of velocity dispersion and surface brightness 
of the JI08 sample. 



3 RESULTS 

In the previous section we have verified the validity of our 
2D decomposition by exploring the 3.6 /xm FP for ellipticals, 
which turns out to be as tight as that observed by JI08 (pri- 
vate communication). Here we study the relations between 
Mbh and the bulge structural parameters listed in Tab. [2] 
and [3]), respectively. To analyze the MeH-bulge scaling rela- 
tions, we adopt three different fitting methods: 

1. a bisector linear regression (Akritas & Ber- 
shady 1996), which uses the bivariate correlated errors 
and intrinsic scatter (BCES) method. Whereas this method 
takes into account the intrinsic scatter, it does not allow 
any determination of it. Hence the intrinsic rms has been 
estimated with a maximum likelihood method assuming nor- 
mally distributed values. 

2. the linear regression FITEXY method as modified 
by T02, that accounts for the intrinsic scatter by adding, in 
quadrature, a constant value to the error of the dependent 
variable in order to obtain a reduced of 1. 

3. a Bayesian approach to linear regression, LIN- 
MIX_ERR (Kelly 2007), which accounts for measurement 



errors, non-detections and intrinsic scatter to compute the 
posterior probability distribution of parameters. 

The MBH-bulge relations we fit are in the following 
form: 

log M. /Mo =a + /3x (x-xo), (2) 

where x is the logarithm of a measured bulge structural 
parameter expressed in solar units and xq is its mean value, 
used to reduce the covariance between the fit parameters. 
As described below, the three methods provide consistent 
results for the MsH-bulge relations. 

Since the LINMIX_ERR is argued to be among the most 
robust regression methods for reliable estimates of the in- 
trinsic dispersion (Kelly 2007), we use it to obtain the final 
results on the Mbh -bulge scaling relations. We still use the 
BCES and FITEXY methods for comparison with previous 
works. 

The fitting results are plotted in Fig. [2] [3] and [5] where 
we present the MeH-^bui, MBH-Mdyn and Mbh-ct relations 
respectively. The Mbh-ct relation is analyzed only to provide 
a consistent comparison with MBH-i^bui and MBH-Mdyn in 
terms of sample and fitting methods. 



M-QYi-hulge scaling relations at 3.6 jim 9 



We exclude from the linear regressions the nine sources 
classified as pseudobulges. This allows us to verify, and pos- 
sibly quantify, whether pseudobulges follow the same scal- 
ing relations as classical bulges. From a visual inspection 
2p 



of Fig, 
10^M( 



pseudobulges with large BH masses (Mbh> 
) follow the relations while those with low BH masses 
deviate significantly from the scaling relations defined by 
classical bulges. Deviant pseudobulges are analyzed in the 
following sections and discussed in Section 4. 

In the following, after describing in details the Mbh- 
bulge relations at 3.6 /xm, we compare them with the pub- 
lished results obtained at shorter wavelengths. Finally in 
Section 3.4 we explore the possible Mbh correlation with 
cr and Re separately as suggested by MH03, Hopkins et 
al. 2007, Graham 2008a. 



3.1 Mbh versus 3.6 /im bulge luminosity 

To construct the MBH-i^3.6,bui correlation, we compute the 
bulge 3.6 /im luminosity using the magnitudes obtained with 
the two-dimensional decomposition (see Section 2.3 for de- 
tails) and reported in Tab. 3. To derive 3.6 /xm luminosi- 
ties in solar units we use the 3.6 /xm Solar absolute mag- 
nitude obtained with a K-band value of 3.3 mag and a K- 
[3.6] =0.05 mag color correction (Allen 1976, see also Bessell 
& Brett 1988, Beh & De Jong 2001). The data for classical 
bulges are fitted with the three linear regression methods 
described at the beginning of this section. We obtain the 
following scaling relation of Mbh with 3.6 /xm bulge lumi- 
nosity (see Tab. 4 for the results from all fitting methods): 



logMBH/M© = (8.19 ±0.06) + 

(0.93=1=0.10) X [log(L3.6,bul/iv3.6,0) 

[rms = 0.38 =1= 0.05] 



11] 



(3) 



All the fit parameters obtained with the three different 
methods are consistent within the errors (see Tab. |4] and 
Fig. [2]). Our results are in excellent agreement with the Mbh- 
^hui,\ obtained at shorter wavelengths. More specifically us- 
ing the BCES and FITEXY results for the comparison, the 
intercept, slope and rms are fully consistent with the K- 
band measurements of MH03, Graham (2007), and H09. We 
find an intrinsic dispersion for the MBH-i^3.6,bui scaling rela- 
tion with rms = 0.35±0.05dex (BCES, FITEXY in Tab|4|, 
consistent with the correlation in the K-band: rms = 0.31 
and 0.36 dex in MH03 and H09 respectively. Optical ob- 
servations in the V-band (Lauer et al. 2007, G09) produce 
slightly higher values for a and (3. Different linear corre- 
lation coefficients can be due to a combination of factors: 
most obviously color corrections, but also incomplete sam- 
ple overlap and the range of Hubble types (G09 fit only early 
types). Furthermore our MBH-^3.6,bui is tighter than the V- 
band correlation whose rms is 0.38 ± 0.09 for early- type 
galaxies only. Moreover extinction can affect the luminosity, 
especially the V-band. 

Figure [2] illustrates that 4 out of 9 pseudobulges 
(NGC3489, NGC3998, NGC4258, NGC4594) are consistent 
with the correlation for classical bulges, while the others 
(CIRCINUS, IC2560, NGC1068, NGC3079, NGC3368) are 
outliers, at more than 4cr below the LINMIX_ERR linear 
regression. The poor statistics for pseudobulges prevent us 



from drawing firm conclusion on the physical nature of this 
behavior. Nevertheless we examine the possible disagree- 
ment also for the other Mbh linear regressions in Section 4. 

3.2 Mbh versus bulge mass 

In this section we first describe two methods to derive the 
bulge mass, and then relate Mbh with Mdyn and M^. Indeed 
one of our aims is to verify whether the 3.6 /xm waveband 
is suitable to investigate the Mbh -bulge mass relations. We 
determine [i) the bulge dynamical mass Mdyn by applying 
the virial theorem (as in MH03, H09), and {ii) the bulge 
stellar mass by calibrating the M/L as a function of the 3.6 
/xm luminosity both with and without V-[3.6] and K-[3.6] 
colors. Thus {i) the bulge Mdyn is computed as follows: 



Mdyn — nRed^ IG, 



(4) 



where G is the gravitational constant, a and Re the veloc- 
ity dispersion (in Tab|2]) and the bulge effective radius (in 
Tab. [3| respectively. It is generally assumed that bulges be- 
have as isothermal spheres and the factor k is usually fixed 
at a value of 3 (Gebhardt et al. 2003) or 8/3 (MH03). How- 
ever, since bulges are well reproduced by Sersic profiles, k 
is in general a function of the Sersic index n. Cappellari et 
al. (2006), using full dynamical models, provide the K(n) 
relation for classical bulges and show that even the use of 
a constant average value of ^ = 5 can provide an accurate 
estimate of Mdyn to within 0.10 dex. We have therefore 
adopted ^ = 5. 

The MBH-Mdyn bulge relation shown in Fig. [3j\ to- 
gether with the three linear regressions. The Mdyn depen- 
dence for the 47 classical bulges is: 



logMBH/Mo = (8.20 =b 0.06) + 

(0.79 =b 0.09) X [log(Mdyn/M0) 
\rms = 0.37 =b 0.06] 



11] 



(5) 



The three linear regressions are consistent within the er- 
rors (See Tab. [4] Fig.jsJV). The MBH-Mdyn relation is rather 
tight with an rms ~ 0.35 dex from LINMIX_ERR, compa- 
rable to the scatter in the Mbh-ct function (see Section 3.3). 
The MBH-Mdyn relations are in excellent agreement with 
the similar K-band analysis in MH03 and H09. The BCES 
regression is also consistent with a sample of 30 galaxies by 
HR04 although their slope is slightly superlinear. 

The position of pseudobulges is the same observed in 
Fig. [2] and described in Section 3.1, where about half of the 
pseudobulges lie more than 4cr below the regression. The re- 
sult for dynamical bulge mass can be compared with those 
for stellar mass to assess wether there are discrepancies be- 
tween the two relations with Mbh- 

To estimate [%%) the bulge stellar mass M^ we have cal- 
ibrated the M/L ratio at 3.6 /xm through color corrections. 
It is known that the M/L ratio is wavelength dependent 
and not constant for a spiral galaxy in any waveband. M/L 
trends have been found to be minimized in the K-band (Bell 
& de Jong 2001, Beh et al. 2003). Published works assume 
a constant M/L with radius, while we are interested in the 
M/L for the bulge alone. The 3.6 /xm wavelength is in prin- 
ciple a good tracer of stellar light for both young and old 
stellar populations as well, thus our data are expected to 
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Figure 2. Scaling relations. The Mbh as a function of 3.6 /im luminosity. The linear regressions are shown as dot dashed blue, dashed 
green and red continuous lines respectively for the BCES, FITEXY and LINMIX.ERR methods and are obtained from classical bulges 
only (47 sources). Pseudobulges are open red squares. 




Figure 3. Scaling relations. The Mbh as a function of bulge dynamical mass (panel a) and bulge stellar mass (panel b). Pseudobulges are 
highlighted as red open squares. In panel (b) the M/L is from equationl6] See Section 3.2 and Fig. [Jlfor details on the M/L calibrations. 
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LINEAR REGRESSION COEFFICIENTS 



log(MBH) - log(x) 


XQ 


Method 


a 




rms 


Ref. 




11.0 


BCES 


8.18 ± 0.06 


0.95 ± 0.10 


0.35 ± 0.02 


- 


^BH-^3.6,bul 


11.0 


FITEXY 


8.21 ± 0.06 


0.91 ± 0.10 


0.35 ±0.05 


- 


11.0 


LINMIX_ERR 


8.19 ± 0.06 


0.93 ± 0.10 


0.38 ± 0.05 


- 




10.9 


BCES 


8.21 ±0.07 


1.13 ±0.12 


0.31 


MH03" 


Lhul,K 


10.9 


FITEXY 


8.29 ±0.08 


0.93 ±0.10 


0.33 


G07 


10.9 


BCES 


8.38 ± 0.05 


0.97 ±0.08 


0.36 


H09 


J-^bul,V 


11.0 


MAX. LIKE.f 


8.95 ± 0.11 


1.11 ± 0.18 


0.38 


G09^ 


10.3 


BCES 


8.41 ±0.11 


1.40 ±0.17 




L07t 




11.0 


BCES 


8.18 ± 0.06 


0.83 ± 0.08 


0.33 ± 0.02 






11.0 


FITEXY 


8.20 ± 0.06 


0.80 ± 0.08 


0.35 ± 0.05 






11.0 


LINMIX_ERR 


8.20 ± 0.06 


79 ± 09 


0.37 ± 0.05 




10.9 


BCES 


8.28 ± 0.06 


0.96 ± 0.07 


0.25 


MH03'^ 




11.0 


BCES 


8.20 ±0.10 


1.12 ±0.06 


0.30 


HR04 




10.9 


BCES 


8.61 ±0.05 


0.88 ±0.06 


0.27 


H09 




11.0 


BCES 


8.15 ±0.05 


0.80 ±0.08 


0.35 ±0.03 






11.0 


FITEXY 


8.17 ±0.06 


0.77 ±0.08 


0.35 ±0.05 






11.0 


LINMIX_ERR 


8.16 ±0.06 


0.79 ±0.08 


0.38 ±0.05 






11.0 


BCES 


8.24 ±0.05 


1.07 ±0.09 


0.32 


H09 




200 


BCES 


8.27 ±0.05 


4.0 ±0.3 


0.32 ±0.04 




Mbh-ct 


200 


FITEXY 


8.30 ±0.05 


4.0 ±0.3 


0.30 ±0.04 






200 


LINMIX_ERR 


8.29 ±0.05 


4.0 ±0.4 


0.33 ±0.04 




Mbh-ct 


200 


FITEXY 


8.13 ±0.06 


4.0 ±0.3 


0.3 


T02 


Mbh-ct 


200 


MAX. LIKE. 


8.23 ±0.08 


4.0 ±0.4 


0.31 ±0.06 


G09^ 



Table 4. Linear regression coefficients of Mbh as a function of the bulge property x, starting from the top of column (1): bulge 3.6 
/im luminosity (equation [sj compared with previous works and with K and V band results; bulge dynamical (equation [sj and stellar 
(equation |8| masses; bulge velocity dispersion (equation [SJ and effective radius (equation |10|. The coefficients are obtained using the 
normalizations in column (2) and with the methods listed in column (3). f G09 use a generalized maximum likelihood method to account 
for upper limits. Columns (4)- (6) list the intercept, slope and rms (in dex) respectively. We compare our results with the works listed 
in column (7). $ Lauer et al. (2007). ^ Coefficients obtained cutting the sample on the basis of the sphere of influence argument for the 
Mbh estimates. ^ Coefficients obtained only for elliptical galaxies. 



trace M/L with a low scatter. Moreover the extent of our 
FP to compact sizes, its tight dispersion valid for all Hubble 
types, and the agreement with the 3.6 /xm FP observed for 
pure elliptical galaxies should guarantee a reliable estimate 
for the M/L of our bulges. 

We assume that Mdyn in bulges is dominated by the 
stellar mass, with a negligible contribution of dark matter 
and gas (Drory, Bender & Hopp 2004; Padmanabhan et al. 
2004). We then consider our sample of spheroids (49 sources, 
see Tab[2]) and fit a linear relation to Mdyn, computed as in 
equation Q, vs iv3.6,bui as: 

log Mdyn = (11.04 ± 0.03) + 

(1.18 ±0.07) X [log(L3.6,bul/i^3.6,0) - 11] 

[rms = 0.13 ±0.04] 

(6) 

This is determined using LINMIX_ERR and is very tight 
with rms ~ 0.10. A tight Mdyn-^3.6,bui relation is not un- 
expected as it is just a different way (i.e. a^Re vs L3.6,bui) 
to express the FP. The Mdyn-^3.6,bui relation is shown in 
Fig.i 

We now include a color correction to test how it af- 
fects the M/L ratio. Only 25 sources are included in this 
exercise, i.e. those whose classical bulges having V- and K- 



magnitudes available in the literature (listed in Tab. 12]). This 
allows a comparison with the above calibration and the fol- 
lowing multi-linear regression: 

log Mdyn = (11.00 ±0.05) + 

(1.35 ± 0.11) X [log(L3.6,bul/i^3.6,0) - 11] 

+ (0.07 ± 0.10) X [{V -K)- 2.54] 
[rms = 0.11 ±0.06] 

(7) 

where we add V-K color correction relative to the sam- 
ple mean V-K color. The fit was performed using the 
multivariate extension of LINMIX_ERR (MLINMIX_ERR, 
Kelly 2007). The zeropoints and slopes of iv3.6,bui in equa- 
tion ^ and ([7f are the same, as well as intrinsic dispersions, 
while the slope of {V — K) is consistent with zero within the 
errors. This, clearly indicates that M/L for bulges is not 
affected by the V-K color and the 3.6 jim. waveband is con- 
firmed to be a better tracer of the stellar light than the 
K-band for both old and young stellar populations. With 
the Mdyn-^3.6,bui relation we have (z) calibrated the M/L 
ratio for the bulge component independently of the Hub- 
ble type of our galaxies, {ii) demonstrated how the 3.6 /xm 
passband traces the stellar bulge mass better than optical 
or near-infrared colors for both young and old stellar popu- 
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Figure 4. M^yn versus L^^.Q^hul- The black line represents the linear regression of equation ([6]). 



lations, and finally {iii) found that V- and K- colors do not 
strongly influence determinations based on L3.6,bui- 
Figure [3)3 shows the Mbh-M^ linear regressions: 

logMBH/M© = (8.16 d= 0.06) + 

(0.79 =b 0.08) X [log(M^/M0) - 11] 
[rms = 0.38 =b 0.05] 

(8) 

where we use equation ([6| to calibrate M^. Equation ([8| is 
just the MBH-i^3.6,bui relation rescaled using equation ^ 
to determine M^, in order to obtain a relation where is 
derived by the host luminosity as done in studies at high 
redshifts (e.g. Merloni et al. 2010). Thus the black hole vs 
bulge mass correlation can be used as the reference scaling 
relation in the local Universe to probe a possible evolution 
of Mbh-M^. This issue is discussed in Section 4. The co- 
efficients and rms of equation (|8| are consistent within the 
errors with the MBH-Mdyn scaling relation (for our choice 
fitting method, see equation ([5|, and Tab. |4|. H09 find a 
steeper slope, but the use of an extinction corrected B-V 
color to calibrate M/L in the K band possibly affect the re- 
sulting Mbh-M^- The M/L calibration presented here pro- 
vide a robust estimate for M^; this issue is pursued further 
in Section 4. 



3.3 A comparison with the Mbh-c" relation 

In the previous sections we have analyzed all the Mbh- 
bulge scaling relations which can be obtained using re- 



sults from our bi-dimensional decomposition. We found {i) 
the same intrinsic dispersion for all the scaling relations 
rms ^ 0.35 dex, and {ii) that pseudobulges have a bi- 
modal distribution with half of them lying on the relations 
within the observed scatter, while those harboring smaller 
BHs (Mbh^ IO^Mq) significantly deviate from the regres- 
sion. Here we compare our Mbh- bulge relations with Mbh- 
cr. A detailed analysis of the Mbh-c" relation is beyond the 
scope of this paper, our aim is just to obtain a consistent 
comparison of Mbh-ct with the other relations in terms of 
samples and fitting method. The classical bulges lie on the 
following correlation: 

log Mbh = (8.29 =b 0.05) + 

(4.0 ±0.4) X log(cr/200 kms~^) 
[rms = 0.33 zb 0.04] 

(9) 

In equation ^ the velocity dispersion is normalized at 
200 km/s according to T02, G09. Figure [s] displays Mbh 
as a function of a and the three linear regressions, with the 
parameters in Tab. [4] 

Although the samples do not fully overlap, our BCES 
and FITEXY results are in good agreement with G09 and 
T02. More interesting is the correspondence, within the er- 
rors, between the MBH-Mdyn (M^, L3.6,bui) and Mbh-c" in- 
trinsic scatters {rms ^ 0.3 — 0.35 dex). Thus at 3.6 /xm is not 
possible to assess whether the bulge dynamics or luminosity 
drive the observed scaling relations. 

As expected, some, but not all, the pseudobulges are 
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Figure 6. Scaling relations. The Mbh is plotted as a function of the bulge effective radius Rq. The linear regressions and symbols are 
as in Fig. |2] 



14 E. Sani et al 



displaced from the inferred regressions. These are CIRCI- 
NUS, IC2560, NGC1068, NGC3079 and NGC3998; there- 
fore there is a good overlap (4/6 where the exceptions are 
NGC3368 and NGC3998) with the deviant pseudobulges 
found in the previous sections. This discrepancy, also found 
in H09 and Nowak et al. 2010, could arise from several dif- 
ferent causes; there could be an additional nuclear compo- 
nent not identified with our image 2D decomposition; or e. 
g. pseudobulges could have a different evolutionary history 
than classical bulges. We defer a detailed discussion to the 
following section. 

3.4 Mbh versus Re, the effective radius at 3.6 /xm 

Mbh has been observed to be separately related to the bulge 
stellar velocity dispersion a and its effective radius Re- This 
has been first observed in MH03 by a partial correlation 
analysis. Figure [6] shows the Mbu-Rs relation obtained with 
our 3.6 fim 2D decomposition, where the LINMIX ERR 
method gives: 

log Mbh = (8.22 zb 0.08) + 

(0.88 =b 0.17) X [\og{Re/kpc) - 0.4] 
[rms = 0.54 zb 0.07] 

(10) 

The three linear regression in Tab. [4] are in good agreement 
within the errors. The intrinsic dispersion is significantly 
larger than in other relations. Nevertheless the possible ex- 
istence of Mbh- Re is intriguing and prompts investigation 
of a possible fundamental plane for BHs, in analogy with 
the galaxies FP. 

It is still unclear whether one MBH-bulge relation is 
significantly tighter than the others. Hopkins et al. (2007) 
explained the lack of a dominant MBH-bulge scaling relation 
because each of them is the projection of the same funda- 
mental plane relating Mbh with two or more bulge proper- 
ties. However, the issue of a BH fundamental plane is still 
matter of debate: Graham 2008a showed that the need for 
a Mbh fundamental plane originates on a displacement of 
barred galaxies with respect to the residuals of the Mbh-ct 
relation. 

If a Mbh FP actually exists, it should be easily observed 
at 3.6 fim. Thus here we follow MH03 and Graham (2008a) 
and compare the residuals of Mbh-ct for our bulges (ex- 
pressed in equation [9| with the effective radius (column 4 in 
Tab. [3]). Figure [t] shows this comparison where we highlight 
the position of barred galaxies and pseudobulges. Hence, we 
check for a possible significant correlation performing both 
the Pearson and Spearman rank correlation test. As shown 
m Tab|5j there is no significant correlation either for the 
entire sample or excluding barred galaxies and/or pseudob- 
ulges. Thus, our analysis at 3.6 /xm does not confirm the 
existence of BH FP. We refer to a forthcoming paper (Mar- 
coni et al. 2011 in prep.) for a detailed study of the origin 
of the FP for supermassive BHs. 



4 DISCUSSION 

The high signal-to-noise images provided by Spitzer /IRAC 
and our improved 2D analysis have allowed us to obtain 



accurate estimates of bulge luminosities and effective radii 
at 3.6 /im, as demonstrated by the tight fundamental plane 
and Mdyn-i^bui relations. We have been able to (i) calibrate 
the M^-Lbui relation, accurate to 0.10 ± 0.05 rms without 
any color correction and (ii) securely define pseudobulges 
through the Sersic index value (n ^ 2) and study their loca- 
tion with respect to the scaling relations for classical bulges. 
We have then obtained the Mbh -galaxy relations for clas- 
sical bulges using our photometry at 3.6 /xm. They are ex- 
pressed in Eqs.[3][5] andjs] and shown in Figs. [2] |3] and [5] 

All the above scaling relations have an intrinsic disper- 
sion of 0.35 dex (from a Bayesian approach to linear regres- 
sion, see Tab[4| and are in perfect agreement with the Mbh- 
ct relation for our galaxies in equation (|9}. Considering the 
methods adopted in literature to perform linear regressions 
that account for the intrinsic scatter (BCES and FITEXY) 
our 3.6 /im correlations are as tight as the Mbh- bulge K- 
band luminosity (MH03, H09) and tighter than in V-band 
(G09). Moreover the MBH-Mdyn and Mbh-M^ agree very 
well with MH03, HR04 and H09, with the new advantage 
of an accurate M/L relation for bulges at 3.6 /xm calibrated 
without any color correction. 

Figure |3] present the Mbh-M^ relation, where M^ is 
obtained by combining Lbui with our calibrated M^-Lbui 
relation. This relation represents the reference Mbh-M^ in 
the local Universe for studies of the MBH-galaxy scaling re- 
lations at higher z, based on M^ derived from luminosity 
measurements (Merloni et al. 2010, Trakhtenbrot & Net- 
zer 2010). The z evolution of MBH-galaxy scaling relations 
has been studied in several papers but remains controversial. 
For instance. Woo et al. (2008) report a significant evolution 
of the BH-bulge relations at z = 0.3 — 0.6, with larger BH 
masses for given bulge cr/luminosity with respect to local 
values. This evolution is however not found by Labita et 
al. (2009) or Shen et al. (2008). Shields et al. (2003) find no 
evolution of Mbh-ct up to 2 but using [OIII] line width 
as a surrogate for a. Indeed, at high redshift a dynamical 
measurement of the bulge mass based on a is difficult due 
to the limited signal-to-noise and spatial resolution. On the 
contrary the total stellar mass of the host galaxy is rela- 
tively easier to estimate from host galaxy luminosities and 
colors (e.g. Peng et al. 2006), or by deconvolving AGN and 
galaxy spectral energy distributions (Merloni et al. 2010) or 
by applying the correlation between the AGN bolometric 
luminosity and Mbh together with the relation between the 
star formation rate and M^ in star forming galaxies (Trakht- 
enbrot & Netzer 2010). The above observational works sug- 
gest that the Mbh / M^ at high redshift is significantly higher 
than the local value. In any case, observations at high red- 
shifts are fundamental to constrain models of BH-galaxy co- 
evolution and to probe different paths followed by galaxies 
during their evolution and the building up of the MBH-bulge 
relations (Lamastra et al. 2010 and references therein). 

In Section 3.2, the tight Mdyn-^3.6,bui relation with neg- 
ligible color corrections justifies the validity of 2D decom- 
position to measure M/L at 3.6 /xm. This (i) unveils the 
capability of our MBH-Mdyn (and thus M^) to be a suit- 
able benchmark to probe the MBH-bulge co-evolution and 
to constrain theoretical models. In fact, in the next decade 
the instruments onboard the James Webb Space Telescope 
(JWST) will increase the possibility to further explore with 
unprecedented accuracy this waveband. 
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Figure 7. Mbh (in Tab[2j residuals about their Mbh-c relation (Eq[9j versus Re. The dashed line represents null residuals. Black 
points are early- and late- type galaxies, blue asterisks are barred galaxies and open red squares are pseudobulges. 



STATISTICAL TESTS FOR Mbh RESIDUALS 



Sample 


Test 


P 


significance 


all galaxies 


Spearman 
Pearson 


0.22 
0.29 


0.10 


no barred 


Spearman 
Pearson 


0.21 
0.29 


0.16 


no pseudobulges 


Spearman 
Pearson 


0.16 
0.24 


0.26 


no pseudobulges, 
no bar 


Spearman 
Pearson 


0.11 
0.20 


0.48 



Table 5. Statistical test for the Mbh residuals about their Mbh-<^ correlation with Re in Fig. |7| Column 1: tests are performed for the 
total sample and sub-samples excluding barred galaxies and pseudobulges. Column 2: the kind of statistic applied to data (Spearman 
rank correlation and Pearson correlation). Column 3: correlation coefficient p. Column 4: and the two-sided significance of p deviation 
from zero in the Spearman rank correlation. The significance is a value in the interval [0.0, 1.0]; a value consistent with zero indicates a 
significant correlation. 



The other {ii) intriguing result found here is related 
to the behavior of pseudobulges. The displacement of pseu- 
dobulges, located below the scaling relations, is reported in 
H09 and Nowak et al. (2010). In §2.4 we followed FDIO 
and identified 9 pseudobulges with Sersic index n ^ 2. Of 
these, CIRCINUS, IC2560, NGC1068, and NGC3398 lie be- 
low all the Mbh scaling relations; NGC3368 deviates from 
MBH-^3.6,bui, -Mdynbut uot from Mbh-c", while the contrary 
happens for NGC3998. The other pseudobulges NGC3489, 



NGC4258, NGC4374, NGC4594 are within the observed 
scatter. As mentioned in §3.3 this discrepancy can arise from 
several causes. We first note that the n ^ 2 is a statistical 
property of pseudobulges; indeed the tails of the Sersic index 
distributions for classical and pseudobulges overlaps around 
n - 2 (see Fig. 9 in Fisher & Drory 2008 and FDIO). Thus, 
not all bulges with a low n would be actual pseudobulges 
from the dynamical point of view (Gadotti 2009). 

Anyway, it is interesting to explore the possible expla- 
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nations once we assume that all the nine galaxies actually 
harbor a pseudobulge. 

(1) Additional nuclear component on small scales 

10 arcsec), such as bars, ovals and lenses are observed not 
only in disks but also in SO galaxies (Laurikainen et al. 2009). 
With the grid method applied for 2D fitting (see §2.3), an 
additional component on small image scales of 10-15 pixels 
can not be taken into account, because (a) it is difficult 
to constrain input parameters while automatically running 
the fits, and (b) the model does not require any additional 
component to match our quality criteria (see §2.3). 

(2) Nowak et al. (2010) define NGC3368 and NGC3489 
as pseudobulges by resolving the nuclear dynamics, and find 
the coexistence of a classical bulge with a pseudobulge, 
where the Mbh seems to better correlate with the classi- 
cal bulge alone. For the same reasons outlined above, it is 
beyond the possibilities of our analysis to deconvolve the 
classical bulge from the pseudobulge component. 

(3) The nine pseudobulges could be in different evo- 
lutionary stages. Montecarlo simulations in Lamastra et 
al. (2010) show how in the buildup of the local Mbu-M^ 
relation, galaxies can follow various evolutionary paths that 
proceed from below the correlation and reach their final po- 
sition passing above it. Thus, the position of each source 
could just be the result of the different evolutionary paths 
being followed. 

(4) Finally, nuclear star clusters (NSCs) may be com- 
mon in those spheroid with 10^ — 10^^ Mq harboring 
a Mbh< 10^ Mq, and may contribute significantly at the 
mass of the central object increasing its value of a factor 
of ^ 2 (Graham & Spitler 2009 and references therein). A 
larger value of Mbh can prevent the coexistence of a NSC, 
because the very central stars should be either disrupted 
by the BH gravity or kicked away by tidal forces. Moreover 
Graham & Spitler (2009) derive a linear regression between 
the nuclear mass (BH plus NSC) and the stellar mass of the 
host spheroid. Thus, the proper quantity to be plotted on 
the y-axis seems to be Mbh plus the NSC mass, at least 
for Mbh< 10'^ Mq. The identification of NSCs requires the 
highest spatial resolutions available with the current instru- 
mentation and can not be reached with Spitzer even for the 
closest of our sources. 

We conclude that the observed displacement of pseu- 
dobulges could arise from observational limitations, true 
evolution, or a combination of the two. 



5 CONCLUSIONS 

In this work we have investigated the scaling relations ob- 
served in the local Universe between Mbh and the structural 
parameters of the host bulges. The analysis is based on a 2D 
decomposition of 3.6 /xm Spitzer/IRAC images of 57 early- 
and late- type galaxies with Mbh measurements. Given the 
well known degeneracy between the bulge Sersic index n 
and effective radius Re, we have adopted a grid-method to 
fit IRAC images and determine n rather than let it freely 
vary. The accuracy of our analysis is verified by the agree- 
ment of our bulges with the mid-infrared fundamental plane 
determination by JI08. Our galaxies extend their 3.6 /im FP 
by doubling the Re range towards smaller radii and none of 
our galaxies deviate significantly from the JI08 FP. This al- 



lows us to reliably (i) calibrate M/L at 3.6 /xm, (ii) study 
the Mbh- scaling relation and (Hi) identify pseudobulges in 
our sample. 

We obtained a tight (rms — 0.10 dex) Mdyn-^3.6,bui 
relation: 

logMdyn/Mo = 11.04 + 1.18 X [log(L3.6,bui/iv3.6,0) - 11], 

which allows us to estimate stellar masses based on luminosi- 
ties at 3.6 /im with high accuracy. The 3.6 /xm luminosity 
appears to be the best tracer of M^ yet found. 

The relations between Mbh, luminosity, masses, and 
effective radius, fitted with a Bayesian approach to linear 
regression are: 

logMBH/M© = 8.19 + 0.93 X [log(L3.6,bul/iv3.6,0) - 11], 

logMBH/Mo = 8.20 + 0.79 x [log(Mdyn/M0) - 11], 

logMBH/M© = 8.16 + 0.79 x [log(M^/M0) - 11], 

logMBH/Mo = 8.22 + 0.9 x [\og{Re/kpc) - 0.4], 

all with the same intrinsic dispersion of rms ^ 0.35 dex 
except for Mbu-Rb which has rms ^ 0.5 dex. Our Mbh-M^ 
relation can be used as the local reference for high redshift 
studies which probe the cosmic evolution of MBH-galaxy 
scaling scaling relations. 

These 3.6 /xm MBH-L3.6,bui, MBH-Mdyn, Mbh-M^ re- 
lations turn out to be as tight as Mbh-c" and, moreover, 
they are consistent with previous determinations from the 
literature at shorter wavelengths. 

We independently identified as pseudobulges those 
galaxies with Sersic index lower than 2 and found 9 sources 
that satisfy this criterion. Of these, 4 pseudobulges lie on 
scaling relations within the observed scatter, while those 
with Mbh lower than 10^ Mq are significantly displaced. 
We discussed the different physical and evolutionary ori- 
gins for such behavior, while considering the presence of 
nuclear morphological components not reproduced by our 
two-dimensional decomposition. 

Finally, we verified the existence of a possible FP for su- 
permassive BHs, relating Mbh with two or more bulge prop- 
erties. We did not find any correlation between the residuals 
of Mbh-ct and the effective radius, showing that our data do 
not require the existence of any BH fundamental plane. 
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APPENDIX A: DETAILS OF THE 2D FITTING 
METHOD 

An important issue for bulge decomposition is related to 
the fitting method. In the case of one dimensional fitting, 
radially averaged surface brightness profiles are fit with 
one or more components. Consequently, information related 
to, e.g., isophotal twists or changes in ellipticity is lost 
and the galaxy profile is poorly defined. Moreover, in one- 
dimensional fitting different morphological components may 
appear to merge smoothly, generating degeneracy and non- 
uniqueness in the profile decomposition. A two-dimensional 
analysis breaks this degeneracy (see, e.g., Peng et al. 2002, 
de Souza, Gadotti & dos Anjos 2004). 

In Section 2.3 we describe in detail our grid method for 
a two-dimensional decomposition of Spitzer /IRAC 3.6 /xm 
images. This is an improved version of the two-dimensional 
method commonly adpted in the literature, where all galaxy 
parameters are left free to vary. In fact, degeneracies also 
affect two dimensional methods: one is related to the func- 
tional form of the bulge profile, "Sersic" profile, in which 
the n and Re parameters are strongly coupled. Moreover, 
the fainter outer parts of galaxies may be affected by er- 
rors on background estimates. In Section 2.3 we also justify 
our choice of the best fit by a quantitative inspection of the 
as a function of the radius. We note that the errors on 
the bulges parameters in Tab. 3 are not one standard de- 
viation errors, rather are given by the discrepancy between 
the best fit values and the ones obtained from the model 
having the closest s, s -\- as or s — as value and the low- 
est possible reduced within Rq^^ according to the criteria 
described in Section 2.3. For each galaxy. Fig. |A1| shows 
the image, 2D model, and residuals. It also shows the trend 
of the source radius for: the best model (i.e. best n 

in red), the fit for the nearest n with low within R^^^ 
maintaining the same sky values (in dark blue, is plotted to 
show the dependence on n at small radii), and the model 
from which we compute the absolute error on free parame- 
ters (in sky blue). Even with the above technique when real 
galaxies are fitted in two dimensions, their surface bright- 
ness in the circumnuclear regions is rarely well reproduced 
only by axisymmetric components. In fact, it is not uncom- 
mon to find inner disks and non-axisymmetric features in 
elliptical galaxies (e.g. Moriondo et al. 1998, Laurikainen et 
al. 2005). Substructure in the circumnuclear regions of both 
elliptical and spiral galaxies (such ovals, disks, lenses) oc- 
curs frequently, and is usually impossible to fit them with 
only axisymmetric components, with or without the addi- 
tion of a bar. This is the origin of nuclear structures in the 
residual images shown in Fig. |A1| F inally, we also note that 



our models (red curves in Fig. Al ) do not require any ad- 
ditional nuclear component to respect the quality criteria of 
Section 2.3, for which, no more than 5% of the nuclear flux 
should be lost. 
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Figure Al. Two-dimensional decomposition. The image, best-fit model and residuals are shown in logarithmic scale. The residuals are 
stretched (±0.25dex) to highlight the finest details. On the right, the x"^ is plotted as a function of the radius (in arcsec). Red: best 
model. Dark blue: the fit with closest n and the lowest within R^l- Light blue: same for the closest sky {s ^(Ts^ see Section 2.3 and 
Appendix A for details). A color version of this figure is available online. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al 



. Continued. 




Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure A2. Continued. 
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Figure A4. Continued. 



30 E. Sani et al 




Figure A5. Continued. 



